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Abstract 

A calculation of the simplest part of the second-order electron self-energy 
(loop after loop irreducible contribution) for hydrogen-like ions with nuclear 
charge numbers 3 < Z < 92 is presented. This serves as a test for the 
more complicated second-order self-energy parts (loop inside loop and crossed 
loop contributions) for heavy one-electron ions. Our results are in strong 
disagreement with recent calculations of Mallampalli and Sapirstein for low 
Z values but are compatible with the two known terms of the analytical Za- 
expansion. 
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The evaluation of the two-loop radiative corrections for hydrogen-like ions with arbitrary 
nuclear charge numbers Z is a challenging theoretical problem important not only for com- 
parison with experimental data for highly charged ions but also for obtaining reliable results 
in the low-Z region. A general review of the present situation can be found in Ref. The 
only part of the two-loop corrections that remains uncalculated up to now is the second- 
order electron self-energy (SESE) part represented by the Feynman graphs in Fig. [l|. The 
first diagram depicted in Fig. |l](a), i.e. the loop after loop irreducible contribution, was 
calculated for Z = 70,80,90,92 |2|] and recently for 1 < Z < 92 ||. This correction is 
separately invariant under any covariant gauge transformation. The corresponding energy 
shift in Ref. |J was called "perturbed orbital" contribution. The results in Refs. agree 
with each other. 

The calculation of the remaining graphs depicted in Figs. |l](b)-(d) is a much more dif- 
ficult problem. All these diagrams have to be evaluated together since only their sum is 
gauge invariant. Part of this sum was calculated recently for Z = 92 with the use of the 
generalization of the potential expansion approach [Q . However, it is not clear whether this 
part is the major one of the total contribution or not. 

Our final goal is also to calculate the remaining SESE corrections given by Figs. [TJ(b) — 
(d) . We are planning to use the renormalization scheme developed in Ref. f| (see also 
Ref. ||) in combination with the partial- wave renormalization approach ]7|J|. Since the 
corresponding numerical calculations are extremely time-consuming (the same holds true 
for the calculations performed in Ref. |4[]), we plan first to apply the approximate approach 
that works only for high Z values and was applied recently for the first-order electron self- 
energy (SE) [PJlOfl. This approach is based on the multiple commutator expansion method 

[n|. 

The original purpose of this work was to test this approximation also on the irreducible 
SESE correction. This means that we had first to calculate the SESE (a) contribution 
exactly and then to compare the result with the approximate one. We had also to elaborate 
the most time-saving procedure compatible with the level of accuracy sufficient for our 
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purpose. Therefore we used the minimal number of grid points and partial waves that could 
guarantee the controlled accuracy within 10%. Working at this level of accuracy and using 
the exact expressions for SESE (a) we found that for the lower Z values our results disagree 
very seriously (more than 50% deviation) with the corresponding results in Ref. ||. The 
comparison will be given below. This disagreement with previous calculations performed by 
Mallampalli and Sapirstein is a crucial point. In Ref. || the breakdown of the perturbation 
expansion in Za was claimed even in the case for Hydrogen. However, perturbation theory 
has been employed intensively in many calculations of radiative effects and thus for the 
purpose of testing QED for weakly bound atomic electrons. Considerable success was made 
in this direction in last years fl2|-|l4||. In contrast to the results and the conclusions drawn in 
Ref. || our numerical results are consistent with the perturbation theory expansion for small 
Z. The fact that we can establish the validity of the Za expansion in low-Z region based 
on exact calculations represents the most important outcome of the present investigation. 
The renomalized SE expression for the atomic state \a) is given by PJT0[] 

AE a = (a\t b (E a )\a) ren = (a\t b (E a )\a) - (a\t f \a) 

oo oo 

= EHsP(£a)|a> ren = £ ((a0 b l \E a )\a) - (a\tf\a)) , (1) 

1=0 1=0 

where Ej and E/ are the bound- and free-electron self-energy operators, and E^ and E^P 
are terms of the corresponding partial-wave expansions. The matrix element of e£ can be 
written as 

(a\E®(E a )\a) = -(21 + 1) £ A na ln |A na | 

x (aWifijii^na^CWn) ■ (n\a%ji(A na r 2 )C l 2 \a) 
- ^(2Z + l)5>gn(£„)A na 

^ n 

x (a77,|a!i l jQ!^(A na r<)nj[(A T i a r>)C' 1 • C l 2 \na), (2) 

where A na = E n — E a , ai^a 2 = 1 — ■ & 2 , ji(z) and rii(z) are the spherical Bessel and 
Neumann functions, respectively, and C\ (i = 1,2) are the standard spherical tensors. We 
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use here the relativistic units h = c = 1 with the fine-structure constant a. The index n 
runs over the whole spectrum of the Dirac equation for the bound electron. 

The matrix elements (a|£? |a) that represent the mass-counterterms in the framework 
of the partial-wave renormalization approach |7]|| read 

(a\tf\a) = ^(2Z + 1) J dp\(p\a)\ 2 J dq A qp In \A qp \ 

x (p|aW«( A gp r i)C'il<z) ■ (q\a%ji(A qp r 2 )C l 2 \p) 
-|(2l + l) J dp\(p\a}\ 2 J dqsgn(E q )A qp 

x (pg|a lM a^ l (A B ,r<)ni(A gp r > )C' 1 • C£|gp), (3) 



where A gp = E q — E p . By means of the symbols |p) and |q) we denote the spherical- wave 
solutions of the free-electron Dirac equation. Integration over p is interpreted as integration 



over the energies E p = ±y77ig + p 2 } where m e is the electron mass and p is the absolute 
value of the electron momentum. The summations over angular quantum numbers are also 
understood. The free-electron wave functions are normalized to 5-function in the energy. 
Equations ([!]) - (|3|) are also valid for arbitrary electronic states \a) and for nondiagonal 
matrix elements of the type (a\Y,(E a )\n) provided that \a) is the ground state. 

The B-spline numerical approach []r5] was used in Refs. P,|T0[ to approximate the sums in 



Eq. (^D and the integrals over q and p in Eq. (||). The number of grid points was N = 140, 
the order of splines k — 9, and the number of partial waves s — 16 ||. The accuracy 



achieved in Ref. compared to the exact Mohrs' results for the point-like nucleus |16 
0.1% for Z = 10 and 0.001% for Z = 92. 



was 



It was observed in Ref. |H| that the terms of Eqs. (|]) - (0) containing ln|A na | are 
dominant for small Z values and that the terms containing sgn(_E n ) are dominant for large 
Z values. From the results obtained in Ref. we can deduce that logarithmic term gives 
99% of the total value for Z = 10 and the sign term yields 95% of the total value for Z = 92. 
The terms analogous to the sign terms in SE can be also found in the expressions for SESE 
corrections. One can try to use the sign approximation for the estimate of the unknown 
parts of SESE for highly charged ions. In this work we will try to test this approximation 



on the SESE (a) correction that can be treated also without any approximations. 
The expression for the SESE (a) correction can be written in the form 

A £irr = g y> (a\^ 1 \E a )\n) rcn (n\t^\E a )\a) ie ^ 
h,h=o 11 E a — E n 

where the summation over n is extended over the whole Dirac spectrum for the bound 

electron and E n = E a term is excluded. For evaluation of the matrix elements in the 

numerator of Eq. @ the formulas @ and (|f) can be applied. Thus in total we have to 

perform 3-fold (or even 4-fold, for counterterms) summations over the spline Dirac spectrum. 

The minimal set of parameters for the numerical spline calculations was chosen to be: 
N = 28, k = 9 and s = 7, while in Ref. N = 50 and 12 < s < 40. However, this minimal 
set in our approach allowed us to keep the controlled accuracy better than 10%. As an 
example, the convergence of our method for the SE correction ([[]) in the ground state for 
Z = 10 is demonstrated in Table |. Both of the partial-wave sequences for odd and even I 
values have been assumed to converge to a common limit. The accuracy of the calculation 
is 7.8% for the minimal basis set. 

We should stress that in our approach, unlike the potential expansion method ||, there 
are no cancellations and no loss of accuracy for small Z values. Still the numerical stability 
becomes poorer in low-Z region. This is a distinct but less dangerous numerical problem. 
The loss of stability for small Z values arises because we employ the spline spectrum gener- 
ated in a large spatial box with the same minimal number of the grid points (N = 28). For 
Z = 1,2, the inaccuracy results to be above the prescribed limit of 10%. 

The results of our calculations of the SESE (a) correction for the ground Is state 
are given in Table [H[ For Z = 70,80,92 values, our results coincide rather well with ones 
in Refs. [@||. The mean deviation is about 1.5%, while the results [0 and || coincide with 
each other within 3 digits. However, for Z = 20 the deviation from Ref. H is about 50% 
and for Z — 10 it is as large as 70%. 

To control the stability of the numerical procedure we compared the results calculated 
with the same iV and s values but with the different order of splines k. In the case of k = 4, 



the deviations from the results obtained in basis set with k = 9 are increased from 1.5% for 
Z = 92 up to 9.5% for Z = 3 and from 30% for Z = 2 up to 50% for Z — 1. According 
to the adopted 10% inaccuracy limit we should consider the results for Z — 1, 2 as unstable 
ones and keep the values only for Z > 3. 

Now we can compare the results for the SESE (a) correction for high Z values given in 
Table [IJ with the results obtained in the sign approximation. The latter arises when we 
retain only the sign terms in all the matrix elements in Eq. The numerical evaluation 
shows that the sign approximation yields 60% of the exact value of SESE (a) for Z = 92. One 
could expect that this approximation may yield results for the other the SESE corrections 
on the same level of accuracy. Such estimates would allow at least to diminish the existing 
uncertainty in the theoretical determination of the Lamb shift for the hydrogen-like uranium 
ions. 

Now let us examine the perturbative nature of QED effects in the \ow-Z region. For 
small Z values we can compare the results of SESE (a) evaluation with the known leading 
terms of Za-expansion of this correction |T3-|I3. We present the result in the standard form 



AET = mJ-) 2 ^GT(Za), (5) 
\7r/ n A 

where n is the principal quantum number of the state a. ^From the Za-expansion calcula- 
tions we know that for small Z values 

G'T(Za) = 2.29953 - ^={Za) \n\Za)- 2 . (6) 



The constant term in Eq. (|6]) was derived in Refs. an d the cubic logarithmic term 



was found in work ||14|| . The results of our calculation of G^J function are given in Table [IT] 
and Fig. ^. In the Fig. ^| the function obtained in Ref. is shown for comparison. 
The results given by the nonrelativistic limit (0) are also plotted. To determine whether our 
results are compatible with the nonrelativistic limit (^[) we tried to use the expression with 
the quadratic logarithmic term 

= 2.29953 - —{Za) \n 3 (Za)~ 2 + C{Za) ln 2 {Za)- 2 . (7) 



To define the coefficient C we use the condition G l { r s = G^* for the different Z values, where 
Gfl is the exact numerical function. For all 3 < Z < 20 we receive nearly the same result 
that after averaging over Z yields C = — 1.0 ± 0.1. The curve corresponding to Eq. (0) with 
the coefficient C = — 1 found from the matching described above is also shown in Fig. |2|. The 
magnitude of the coefficient reveals that our results are consistent with the Za-expansion 
perturbation theory |l2|-p^]. The more detailed comparison with the results obtained in 
Ref. |§ shows that the total difference comes from the sum over negative states in Eq. (f|) 



T7fl . The reason for this discrepancy could be not only the difference in the spline spectrum 



but also the difference in the evaluation of the off-diagonal self-energy for the negative states. 
The latter is very difficult to compare explicitly since the methods used for this evaluation 
in this work and Ref. || are quite different. 
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FIGURES 




(a) (b) (c) (d) 

FIG. 1. The second-order electron self-energy Feynman graphs. The double solid line denotes 
the electron in the field of the nucleus, the wavy line denotes the photon. The double line with the 
bar denotes the electron propagator with the reference state excluded from the summation over the 
Dirac spectrum. The symbol in (b) denotes the derivative of SE graph over an energy parameter 
E in the bound electron propagator. The graphs (a) and (b) correspond to the irreducible and 
reducible parts of the loop after loop contribution, the graph (c) corresponds to the loop inside 
loop contribution, and the graph (d) corresponds to the crossed loopes contribution. 
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FIG. 2. The function G"(Za) in the different calculations. 
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TABLES 

TABLE I. Convergence of the partial J2i Ai^il contributions in the first order self-energy AEi, 
evaluation for Z = 10. The number of the grid points was chosen to be N = 28. 

even odd AE ls (eV) 

1 Ej'=o AE u ( eV ) 1 E|'=i AE u ( eV ) This work Ref - [0 

0.1269 1 0.1769 0.1688 0.1566 

2 0.1615 3 0.1715 
4 0.1659 5 0.1702 
6 0.1672 



12 



TABLE II. Comparison of the SESE (a) correction for the ls-ground state with previous cal- 
culations. 



AJ3& (eV) GT S 



z 


Ref. H 


Ref. § 




This work 




Ref. i 


This work 


3 




-0.6237 x 10" 


-7 


-0.2913 x 10- 


-7 


-4.50 


-2.101 


4 




-0.2786 x 10" 


-6 


-0.1351 x 10- 


-6 


-4.77 


-2.311 


5 




-0.8792 x 10- 


-6 


-0.4431 x 10- 


-6 


-4.931 


-2.485 


6 








-0.1153 x 10- 


-5 




-2.599 


7 




-0.4808 x 10- 


-5 


-0.2584 x 10- 


-5 


-5.016 


-2.694 


8 








-0.4972 x 10- 


-5 




-2.659 


9 








-0.8903 x 10- 


-5 




-2.642 


10 




-0.2796 x 10- 


-1 


-0.1483 x 10- 


-4 


-4.9016 


-2.601 


20 




-0.7525 x 10- 


-3 


-0.4688 x 10" 


-3 


-4.1217 


-2.568 


30 








-0.3454 x 10- 


-2 




-2.491 


50 








-0.4407 x 10- 


-1 




-2.472 


70 


-0.2283 


-0.2282 




-0.2314 




-2.3804 


-2.413 


80 


-0.4474 


-0.4472 




-0.4512 




-2.3923 


-2.413 


92 


-0.9712 


-0.9706 




-0.9599 




-2.581 


-2.553 
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